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We consider the possibility of energy being exchanged between the scalar and matter fields in 
scalar-tensor theories of gravity. Such an exchange provides a new mechanism which can drive 
variations in the gravitational 'constant' G. We find exact solutions for the evolution of spatially 
fiat Friedman-Robertson- Walker cosmologies in this scenario and discuss their behaviour at both 
early and late times. We also consider the physical consequences and observational constraints on 
these models. 
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CN ■ I. INTRODUCTION 

^ ' Scalar-tensor theories of gravity provide a convenient framework within which to model space-time variations of the 
^> , Newtonian gravitational constant, G. They feature a scalar field, (j>, which is non-minimally coupled to the space-time 
curvature in the gravitational action. It is this scalar, or more usually its reciprocal, that drives variations in G. 
, Non-minimally coupled scalar fields arise in a variety of different theories, including Kaluza-Klein theory jlj], string 
theories (3| and brane- worlds The same mechanism that creates a scalar field non-minimally coupled to the 
' curvature in these theories can also lead to a coupling between the scalar and matter fields. This coupling manifests 
^ itself through the matter Lagrangian becoming a function of 4>. The possibility of such a coupling is usually neglected 
^sO ill the literature, where the matter Lagrangian is a priori assumed to be independent of (j). It is the possibility of a 
T-H , coupling between the scalar and matter fields in scalar-tensor theories that will be the subject of this work. 

' The introduction of a coupling between cj) and matter greatly enlarges the phenomenology of the theory. Potentially, 
! this allows greater variability of G in the early universe whilst still satisfying the solar system bounds on time-varying 
■ G p. We will consider spatially- flat Friedmann-Robertson- Walker (FRW) cosmologies and investigate the extent to 
I which G can vary when energy is exchanged between the (j> field and ordinary matter. As well as giving a window into 
the four-dimensional cosmologies associated with higher-dimensional theories, we hope that this direction of study 
might also be useful in understanding why the present value of G is so small compared to the proton mass scale 



X . {Grrip^ ~ 10 ^^). The direct exchange of energy between (j> and the matter fields offers a non-adiabatic mechanism for 



X 



G to 'decay' towards its present value from a potentially different initial value. There have been a variety of studies 
which investigate the drain of energy from ordered motion by entropy generation, due to bulk viscosity |5j or direct 
decay 0,0] or energy exchange Q , but few studies of the drain of energy by non-adiabatic processes from a scalar field 
that defines the strength of gravity [^ ITollTTI| . This creates a range of new behaviours in scalar-tensor cosmologies. 
In considering a coupling between (f> and matter we are forced to reconsider the equivalence principle. The energy- 
' momentum tensor of perfect-fiuid matter fields will no longer be covariantly conserved and the trajectories of test- 
particles will no longer follow exact geodesies of the metric. These violations of the experimentally well verified weak 
equivalence principle exclude most possible couplings between (j) and matter 0. Such violations are not necessarily 
fatal though. We show that whilst energy-momentum is not separately conserved by the matter fields there is still an 
exact concept of energy-momentum conservation when the energy density of the scalar field is included. Furthermore, 
the non-geodesic motion of test particles is only problematic if the coupling increases above experimentally acceptable 
levels as the Universe ages. The theory we consider is still a geometric one and it remains true that at any point on 
the space-time manifold it is possible to choose normal coordinates so that it looks locally flat, ensuring that it is 
always possible to transform to a freely- falling frame in which the effects of gravity are negligible (up to tidal forces). 
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II. FIELD EQUATIONS 

The simplest scalar-tensor theory is the Brans-Dicke theory (l3|, defined by the Lagrangian density 

£^0i?_ ^9-09^0 +167r£„[g,b;*] (1) 

where lo is the Brans-Dicke coupling constant, i? is the scalar curvature of space-time, and £m is the Lagrangian 
density of the matter fields, denoted by As ^ oo this theory reduces to general relativity and G ~ becomes 
constant (for exceptions see It is an important feature of these theories that Cm is independent of </>. This 

ensures that the matter fields do not interact with the scalar field directly and therefore that the energy-momentum 
tensor, T"-^ , derived from Lm is conserved {T°-^.f^ = 0). 

This conservation of T"*", whilst appealing, is not absolutely necessary in deriving a theory in which G can vary. 
There are numerous examples where one might expect r"**,^ ^ 0. For example, when considering two fluids the 
energy-momentum tensor of each fluid is not separately conserved unless the fluids are completely non-interacting. It 
is only required that the energy-momentum being lost by one of the fluids is equal to the energy-momentum being 
gained by the other. 

In what follows we will consider the scalar field and the matter fields as two fluids (or more than two fluids if there 
is more than one matter fluid present) and introduce a transfer of energy and momentum between them. Such an 
interaction can be introduced by allowing Cm to be a function of (j) and will change the nature of the resulting FRW 
cosmologies. 

The field equations are derived from by extremizing the corresponding action with respect to the metric. Defining 
T"'' for the matter in the usual way, the field equations take their standard Brans-Dicke form independent of the 
presence of interactions between and matter, and the Einstein tensor is given by 

G^" - ^(0;V;' - ^5'^V;c0;^) + ^(</>;'"' " S'^'D^) + jT,f . (2) 

The scalar-field propagation equation and matter energy-momentum conservation equations are 

^ {2uj + 3) i2u; + 3) (j)," ^ ' 

T^^.b = cr" (4) 

where T is the trace of the energy-momentum tensor and cr" is an arbitrary vector function of the space-time coor- 
dinates cc'' that determines the rate of transfer of energy and momentum between the scalar field (p and the ordinary 
matter fields. The precise form of cr" depends on the detailed form of the interaction between the scalar and matter 
fields in Cm- For example, a conformal transformation of the form gab A^{4>)9ab from a frame in which T°'\ = 
gives 

a TdA 

This particular choice of energy transfer can be interpreted as a space-time variation of the rest masses of matter 
described by Cm- For the moment, we consider the case of more general interactions by leaving cr° as an arbitrary 
function. Later, we will consider specific forms of <t° that allow direct integration of the field equations. 

We specialise the metric to the spatially flat, isotropic and homogeneous FRW line-element with expansion scale 
factor a{t): 

ds^ ^~dt^ + (t) {dr^ + {d0^ + sin^ Odc/)^ )). (5) 
Substituting this metric into the field equations lO, © and 10} gives the generalised Friedmann equations: 

(6) 

(M=3«((2-,V + l) + 3^^.o (,) 

^=«(4-3,)-2^^.° (8) 
pV P (j) 

P + l^P^^"" (9) 
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where we have defined a comoving volume V = a"^ and a constant a = j^^^', the energy-momentum tensor is 
assumed to be a perfect barotropic fluid with density p and pressure p which are linked by a linear equation of state 
p — {"f — and over-dots denote differentiation with respect to the comoving proper time, t. It is this set of 
differential equations that we need to solve in order to determine the evolution of a(t) and G oc (f>{t)~^ in cosmological 
models of this type. 



III. TRANSFER OF ENERGY AND ENTROPY 



The conservation of energy and momentum as well as the second law of thermodynamics are of basic importance 
to physics. In considering an interaction between a gravitational scalar field and the matter fields 'J it is, therefore, 
necessary to investigate the extent to which we can consider energy and momentum to be conserved and the second 
law to be obeyed. 

When we consider the thermodynamics of an exchange of energy between the scalar field and matter it is useful to 
define an effective energy density, p^, for the scalar field 0. Defining 

(10) 



the scalar-field propagation equation ((TJ can then be rewritten as 

R ■ 



(11) 



Comparison of this equation with Q shows that (j) acts as a fluid with equation of state ^ — 2 (p^ — p^). The 
two terms on the right hand side of this equation act as sources for the energy density p^. The first is the standard 
Brans-Dicke source term for the scalar field and the second, cr'^{t), is new and describes the energy exchange between 
(j) and the matter fields. It can be seen that the second term is exactly the opposite of the source term in equation 
and it is in this sense that the total energy is conserved in this theory. 
It is also useful to consider the entropy. Contracting the divergence of the energy-momentum tensor with the 
comoving four-velocity C/" we obtain 

= U''P;a + U''{{p+p)UaU''),b 

= uy,a-{{p+p)u\b, 

where, in the last line, we have used the normalisation W^Ua — —1 and {U°'Ua);b = U°'-bUa + U'^Ua-b = 0. Defining 
the particle current by N°' = nU"", where n is the number density in a comoving Lorentz frame, this expression can 
be rewritten as 



P;a - n 



(P+P) 



iP+P) 



P - 

n 



where we have used the conservation of particle number, N'^-a — 0. Recalling the first law of thermodynamics. 



edS = pdV + dE^pd[- 



where Q is the temperature and S is the entropy, we now get 
or, making use of our assumption of spatial homogeneity, 



di^ 
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This tells us that as energy is transferred from (f> to the matter fields the entropy of the matter fields increases, as 
expected. Conversely, the matter fields can decrease their entropy by transferring energy into (j). 

Unfortunately, there is currently no known way of defining the entropy of a non-static gravitational field so it is 
not possible to perform an explicit calculation of the entropy changes in (j) and gab- We can only assume that if the 
Universe can be treated as a closed system, and the exchange of energy is an equilibrium process, then the entropy 
that is lost or gained by the matter through this exchange will be gained or lost by the gravitational fields. This 
direct interaction of the matter with then allows an additional mechanism for increasing or decreasing the entropy 
of the matter content of the Universe. 



IV. GENERAL SOLUTIONS 



It is convenient to define a new time coordinate r by 



dr = pVdt (12) 

and to re-parametrise the arbitrary function a" by 



- p^V^X' (13) 



where a prime denotes differentiation with respect to r and A(r) is a new arbitrary function. This re-parameterisation 
of the interaction is chosen to enable a direct integration of the field equations and does not imply any loss of generality, 
as A is an arbitrary function. The field equations {Tj) and ^ can now be integrated to 

p(j)VV' = 3q((2 - 7)cj + 1)(t - Ti) + 3aA (14) 
pV'^(j)' ^a{'i-i-l)T + aT2-2a\ (15) 

where ri and T2 are constants of integration. We have a freedom in where we define the origin of r and can, therefore, 
absorb the constant ri into r and the definition of T2 by the transformations r — > r -I- ti and r2 ^ T2 — (4 — 37)ti. 
It can now be seen from (I15|l that (f)' is sourced by three terms. The first corresponds to the source term in and 
can be seen to disappear for 7 = 4/3, as expected for black-body radiation. The second term is constant and is 
the contribution of the free scalar to the evolution of (j)] it is this term which distinguishes the general spatially-fiat 
Brans-Dicke FRW solutions j3| from the power-law late-time attractor solutions The third term is new and 

gives the effect of the energy transfer on the evolution of cf). This term is dependent on the arbitrary function A, which 
specifies the interaction between cj) and the matter fields. 

The problem is now reduced to solving the coupled set of first-order ordinary differential equations H14|) and \ib\ 
with the constraint equation (jSJ. The remaining equation ((HJ is rewritten in terms of r and A as 

p V (j) 

and can be solved for p once V and cf) have been found for some A. 

We can decouple the set of equations p4(l and H15|l by differentiating p5|l and substituting for l|14|l to get the 
second-order ordinary differential equation 

((4 - 37)t + r2 - 2A) ^ - ((4 - 37) - 2A') ^ 

= - [((4 - 37)t + T2- 2A)A' + 3(2 - 7)(((2 - 7)c^ + l)r + A)] 
which can be integrated to 



where 



£ _ (4 - 37)t + T2-2\ 
(j) ~ {At'^ +Bt + C) ' 

A ^ 3j^Lj/2 - 37(1 + 2uj) + (5 + 6cj) 
B = T2 + {A- 37)A 
C = -A^ -t- At2 -I- D 



(17) 
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and D is a constant of integration. The three source terms for 0' appear in the numerator on the right-hand side of 
equation H17|) . The equations H14|) and (|15|l can now be combined to give a! ja in terms of (\>' l(\> as 

o! ^ ((2-7V + l)r + A 

a (Ar2 + Sr + C) ^ ' 

where A, B, C and D are defined as before. The constant D can be set using the constraint equation Using H14|l 
and H18|l we obtain the expression 

pV^(j) = a{AT'^ +Bt + C). 

This can then be substituted into ^ which, in terms of r and a, gives the generahsed Friedmann equation: 

a'Y a' (j)' uj (3 + 2w) 



3 - +3 



a J a (p 2 J [At"^ + Bt + C) 
Substituting (|17|l and (|18|) into this we find that 

D 



2(3 + 2w)" 

V. PARTICULAR SOLUTIONS 

If A is specified in terms of r we now have a set of two decoupled first-order ordinary differential equations for the 
two variables a and 0. It is the solution of these equations, for specific choices of A(t), that wc give in this section. 

A. A(r) = ci -I- C2T 

A simple form for A that allows direct integration of equations and (|18|l is the linear function A = ci +C2T. From 
equations H14|l and H15(l it can be seen that the constant ci can be absorbed into ri and T2 by simple redefinitions. 
The equations (flTI and fH^ then become 

_ (4 - 37 - 2c2)t -f- T2 
^ {At^ + Bt + C) 

_ ((2 - 7)^ + 1 + C2)r 
a ~ (At^ + Bt + C) 

where A = A—c\, B = B+C2T2 and C — D. The solutions of these equations depend upon the sign of the discriminant 

A = - AAC. (19) 
For the case A = 0, there exist simple exact power-law solutions 

2(2-7)^ + 2 + 2c2 



a(T) oc r^"'""-«"'(i+^-'+=<^+«"'-2'=2 (20) 

2(4-37)-4c2 

(/)(t) oc '5^<i+=-'+=^<^+«"'-2'=2 . (21) 



Substituting these power-law solutions into (|16|l we can obtain the corresponding power-law form for p 

4e2(2-37-c2)-67(l + (2--,)i^) 
p ^ ^372„_67(l + 2i..)+2(5+6u)-2c| ^ 

The relationship between r and the cosmological time t can now be obtained by integrating the definition dr — pa^dt 
given in eq. (|12|l . This gives l|20|l and H21|l in terms of t time as 

2 + 2(2-7)^ + 2c2 

a{t) ~ i4 + 37c^(2-7)-2c2(7-67-c2) (22) 

2(4-37)-4c2 

(j){t) ~ t4 + 37c.(2-7)-2c2(7-67-c2) , (23) 
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The condition required for the occurrence of power-law inflation is obtained by requiring the power of time in equation 
(|23|l to exceed unity (for the case with out energy transfer see refs. For uj > —3/2 we always have A ^ 0, 

and the case A > possesses the exact solutions 



a(r) = aoiAr^ + Bt + C) ? 



-B((2--f)i>, + l + c2) 



2 At + ]3-^fK^ 



{t) = MAt' + Bt + C) — 



(4-3t-2c2)B-t2^ 



2 At + B-yfK , 



(24) 



(25) 



where ag and 0o are constants of integration. For lo < —3/2 we have A ^ 0, and the case A < has the exact 
solutions 



i(t) = ao(AT^ + Sr + C) exp 



((2-7)w + l + C2)B , _i/'B + 2Ar 

— L&.I1 I 



A^ 



(r) = 00 (ir^ + Br + C*)* '2^^"' exp 



2r2i - 2(4 - 37 - 2c2)B . ( B^ 2At\ 

tan I . 



AJ^ 



-A j 



(26) 
(27) 



These solutions have the same functional form as those found by Gurevich, Finkelstein and Ruban for Brans-Dicke 
theory, in the absence of energy exchange (A — constant), and reduce to them in the limit C2 0. The behaviour of 
these solutions at early and late times will be discussed in the next section. 

B. A(r) = car", n / 1 

We now consider forms of A(r) that are more general than a simple linear function of r. Making the choice A — c^t" , 
where n 1 and C3 is constant, and setting the free scalar component to zero (r2 = 0), we find that H17|l and H18|l 
can be integrated exactly. The form of the solutions again depends upon the roots of the denominator. For real roots 
we require u > —3/2, for which we find the solutions 



a(r) = uqt' 



±2c3r"-i T (4 - 37) ± (2 - 7)^3(3 + 2c^) 



3 + 3(2-7)^-^3(3 + 20.) 
3k(,.-1) 



±2c3r"-i T (4 - 37) T (2 - j)^/3{3 + 2u;) 



3 + 3(2 — ,)lo + ^3(3 + 2lj) 
3k(,i-1) 



(T) = <PoT 



2(4-3t) 



±2c3r"-i T (4 - 37) ± (2 - 7)^3(3 -I- 2c^) 



(4-3t) + (2-7)V3(3 + 2i^) 



±2c3r"-i T (4 - 37) T (2 - 7)^3(3 + 2c^) 



(4-3t)-(2-t)V3(3 + 2i^) 



(28) 



(29) 



where k = 2(5 — 37) + 3(2 — For a denominator with imaginary roots we require uj < —3/2, for which we find 

l+(2-7)^ 



(r) = ao ±2(4 - 37)c3ri-" ± Kr^^i"") =f 2c] 



xexpi^^gi±Mtan-/ (4 - 37) - 2c3t 



3At(l-n) \^(2-7)^-3(3 + 2w)^ 

(4-37) 



(r) = 00 ±2(4 - 37)c3ri-" ± Kr^^^-") ^ 2c^ 



j 2(2-7)V-3(3 + 2c^) ^ _J (4 - 37) - 2c3<'^ 
X exp < f: ^ tan ' 



k(1 — n) 



(2- 7)^-3(3 + 2l^) 



(30) 



(31) 



with K defined as above. The ± and =p signs here indicate that there are multiple solutions that satisfy the field 
equations. These signs should be chosen consistently within each set of square brackets (solutions H28|) and H29|l can 
have upper or lower branches chosen independently in each set of square brackets, as long as a consistent branch is 
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taken within each separate set of square brackets). The physical branch should be chosen as the one for which the 
quantity in brackets remains positive as r — *■ oo, so ensuring the existence of a positive real root in this limit, 

These solutions display interesting new behaviours at both early and late times, which will be discussed in the next 
section. 



VI. BEHAVIOUR OF SOLUTIONS 



These exact solutions for a(T) and ff>{T) can now be analysed at early and late times. 



A. A(r) — ci + C2T 



At late times, as r ^ oo, the solutions H24() - (|27() all approach the exact power-law solutions (|20|) an d (|21l) . It can be 
seen that these solutions reduce to the usual spatially flat FRW Brans-Dicke power-law solutions [l^ in the limit that 
the rate of energy transfer goes to zero, C2 — *■ 0. It can also be seen that these solutions reduce to the spatially-flat 
FRW general relativistic solutions in the limit lo oo, irrespective of any (finite) amount of energy transfer. 

The early-time behaviour of these solutions approaches that of the general Brans-Dicke solutions without 
energy transfer. Generally, we expect an early period of free-scalar-field domination except in the case T2 = 0, in 
which case the power-law solutions (|20f) and H21I) are valid right up to the initial singularity. For uj > —3/2, the 
scalar-field dominated phase causes an early period of power-law inflation. In this case there is always an initial 
singularity and the value of the scalar field diverges to infinity or zero as it is approached, depending on the sign of 
T2. For UJ < —3/2 there is a 'bounce' and the scale factor has a minimum non-zero value. In these universes there is 
a phase of contraction followed by a phase of expansion, with no singularity separating them. Solutions of this type 
were the focus of (2j where the evolution of through the bounce was used to model the variation of various physical 
constants in such situations. The energy exchange term does not play a significant role at early times in these models. 
The asymptotic solutions as the singularity (or bounce) is approached are the same as if the energy exchange term 
had been neglected, and are given by , up to the absorption of ci into T2 previously described. 



B. 



A(r) = car", n / 1 



The behaviour of the solutions H28|l - H31|l depends upon the signs of n— 1 and C3, as well as on the sign of w-l-3/2. For 
illustrative purposes we will consider the radiation case 7 = 4/3 which is appropriate for realistic universes dominated 
by asymptotically-free interactions at early times. 

For n > 1 it can be seen that the late-time attractors of solutions (|28f) -(l31 |l are a constant and (p — * constant as 
T 00, for both Lu > —3/2 and uj < —3/2. At late times these universes are asymptotically static; the evolution of 
the scale-factor ceases as r — > 00 and both (j) and p become constant. Further analysis is required to establish whether 
these static universes are stable or not (we expect them to be stable as no tuning of parameters or initial conditions 
has been performed to obtain these solutions). 

The early-time behaviour of solutions with n > 1 depends upon the sign of C3 as well as whether uj is greater or less 



than —3/2. We will consider first the case of lo > 



see that a ^ 00 as r*^ 



Tq (where Tq = ((2 



3/2. For C3 > 0, we see that a 
7)^3(3 



2lu) + (4 - 



00 as T 
37))/2c3 and Tq = 



Tq ; for C3 < we 



-((2- 7)^3(3 + 2l^)- 

(4 — 37))/2c3). For n > 1 and ui > —3/2 we therefore have the generic behaviour that a ^ 00 at early times and a 
constant at late times. The behaviour of a at intermediate times varies in form depending on the sign of C3, as can 
be seen in Figures 1(a) and 1(b) The asymptotic form of for n > 1 and uj > —3/2 depends critically on the sign of 



Tg , whereas for C3 < it can be seen that 



C3. For C3 > it can be seen that ^ as r" 

The behaviour of in these two cases is illustrated in Figures 1(c) and 1(d) 



00 as T 



We now consider the early-time behaviour of solutions with n > 1 and uj < —3/2. It can be seen that a — > as 
T — > irrespective of the sign of C3, so that we find the generic behaviour a — s- at early times and a — > constant at 
late times (this is in contrast to the standard theory where an initial singularity is avoided when uj < —3/2 ). Again 



the behaviour of a at intermediate times is dependent on the sign of C3, as can be seen from Figures 2(a) and 2(b) 
As T — > we see that (f> has a finite non-zero value and is either increasing or decreasing depending on the sign of C3 
This behaviour is shown in Figures 2(c) and 2(d) 




(c) C3 = 10 (d) C3 = -10 

FIG. 1: The time evolution of a and for n = 2, w = 10 and 7 = 4/3. 




(c) C3 = 10 (d) C3 = -10 

FIG. 2: The time evolution of a and </> for n = 2, a> = — 10 and 7 = 4/3. 



(a) C3 = 10 



(b) C3 = -10 




(c) C3 = 10 



(d) C3 = -10 



FIG. 3: The time evolution of a and for n = 0, w = 10 and 7 = 4/3. 




(c) C3 = 10 (d) C3 = -10 

FIG. 4: The time evolution of a and </> for n = 0, w = — 10 and 7 = 4/3. 
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It remains to investigate the nature of the solutions with n < 1. At late times we see that 



2 + 2(2 — 
a T-2(5-3T) + 3(2-T)i:c 

2(4-3t) 
6 t-2(5-3t) + 3(2-t)^i. 



as r ^ oo, irrespective of the sign of C3 or the value of cu. These late-time attractors are the flat FRW power-law 
Brans-Dicke solutions [l9l | which reduce to the standard general-relativistic solutions in the limit w — > 00. 

The early-time behaviour when n < 1 depends on the sign of C3 and the sign of + 3/2. We consider first the case 
u > —3/2. In this case it can be seen that a ^ as t"~^ — > Tq or r" ~^ ~lUh.^ ^'-'^ C3 > or C3 < respectively. 
The behaviour of a for both C3 > and C3 < is shown in Figure p(a)| and 3(b) The behaviour of (j) at early times 
depend s on t he s ign o f C3 and goes to 00 for C3 > or to for C3 < 0. The behaviour of 4> in these cases is shown in 
Figures 3(c) and 3(d) For uj < —3/2 the scale factor a contracts to a finite, but non-zero, minimum value and then 
expands. The exact form of the minimum depends on the values of n, ui and C3, but it is interesting to note that odd 
values of n produce symmetric bounces and even values of n produce asymmetric bounces, as illustrated in Figures 
4(a) and 4(b) The evol ution of c j) thro ugh these bounces is smooth with a time direction prescribed by the value of 
C3, as shown in Figures 4(c) and 4(d) (increasing for C3 > and decreasing for C3 < 0). The effect of changing the 



sign of C3 is seen to be a mirroring of the evolution of a and (p in the y— axis. 



VII. PHYSICAL CONSEQUENCES 

The solutions found in the previous sections are of physical interest for a number of reasons. The transfer of 
energy and momentum between a non-minimally coupled scalar field (j) and matter fields is a prediction of a number 
of fundamental theories of current interest, including string theories, Kaluza-Klein theories and brane-worlds. The 
cosmologies produced by such an interaction, therefore, should be of direct interest in the consideration of these 
theories. Furthermore, the solutions we have found display modified behaviour at both early and late times. The 
investigation of modified theories of gravity at early times is of particular interest as it is in the high-energy limit that 
deviations from general relativity are usually expected. Modified behaviour at late times is also of interest as it is at 
these times that we can make direct observations which can be used to constrain deviations from the standard general- 
relativistic model. We will now summarise the behaviour of the solutions found in the previous section, highlighting 
the physically significant results and constraints that can be placed on the theory by observations. 

For the case of A linear in r it was shown that the late-time attractors of the general solutions are no longer the 
power-law solutions of the Brans-Dicke theory, but are given by equations (|22|l and (|23|l . These attractors are of 
special interest as they have a simple power-law form that reduces to the general relativistic result in the limit uj 00 
and to the Brans-Dicke result in the limit C2 0. Observations of cosmic microwave background anisotropics and 
the products of primordial nucleosynthesis will therefore be able to constraint any potential late-time deviations of 
this kind, and hence the underlying model. The process of primordial nucleosynthesis in scalar-tensor theories has 
been used by a number of authors to place constraints on the coupling parameter uj{(f>) |2lll23 . 12M 12^ . I25L l2^ l27j | . In 
these studies the different value of G during nucleosynthesis causes the weak interactions to freeze out at a different 
time and hence the proton to neutron ratio at this time is different to the standard case. This modification causes 
different abundances of the light elements to be produced, which can be compared with observations to constrain 
the underlying theory. Studies of this kind usually assume G to be constant during nucleosynthesis, which will not 
be the case when e nerg y is allowed to be exchanged between (j) and the matter fields. The effects of a non-constant 
G were studied in |27l |. A similar study would be required to place constraints upon the parameters C2 and uj in 
this theory. The cosmic microwave background power spectrum has also often been used to constrain scalar-tensor 
theories of gravity [2^ [2^ [s^, |^ IH, HI] . In these studies the redshift of matter-radiation equality is different from its 
usual general relativistic value due to the modified late-time evolution of the Universe. This change in the redshift of 
equality is imprinted on the spectrum of perturbations as it is only after equality that sub-horizon scale perturbations 
are allowed to grow. The main effect is seen as a shift in the first peak of the power-spectrum, which can be compared 
with observations to constrain the theory. Again, the late time evolution of the Universe is modified from the usual 
Brans-Dicke case by the energy exchange that we consider, so that the previous constraints are not directly applicable. 

For the case of a non-linear power-law exchange of energy, described by A cx r" , the late-time evolution of a and 
(j) can be significantly modified. For n > 1, the solutions do not continue to expand eternally, but are attracted 
towards a static state where the time-evolutions of a, (j) and p cease. For n < 1 the generic late-time attractor is 
the power-law solution of a flat FRW Brans-Dicke universe. It appears that theories of energy exchange with n > 1 
are ruled out immediately by observations of an expanding universe whilst the case of n < 1 is subject to the same 
late-time constraints as the standard Brans-Dicke theory |m I2I I2I I2I I2I I2I I2I I2I I3I IH l3^ 
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FIG. 5: The time evolution of a and (j) for n = —2, uj = —10, ca = 10 and 7 = 4/3. 

It remains to investigate the physical consequences of the early-time behaviour of our solutions. For A linear in r 
these solutions approach those of the standard Brans-Dickc theory as either the initial singularity or the minimum 
of the bounce are approached, according to the sign of w — 3/2. The physical significance of this behaviour has been 
discussed many times before, usually focusing on the avoidance of the initial singularity and the inflation that can 
result from the presence of the free component of the scalar field. 

For A c>c t" the early-time behaviour can be significantly changed from that of the standard theory. For n > 1 the 
scale factor a either approaches infinity or zero, depending on our choice of u and C3, as previously described. For 
the more realistic case of n < 1 the evolution of a at early times either undergoes a period of rapid expansion or a 
non-singular bounce, depending on whether to is greater or less than —3/2. This behaviour is similar to that of the 
general solutions of the standard theory, but in this case the free scalar-field-dominated epoch has not been invoked 
and there is more freedom as to the exact form of the evolution. For example, with a suitable choice of parameters 
it is possible to create a universe that contracts and then is briefly static before 'bouncing' and continuing on to its 
late-time power-law evolution. This is shown in Figure|31for the case uj = —10, n — —2 and C3 = 10. (It is interesting 
to note that Peter and Pinto-Neto remark that a static period followed by a bounce could potentially produce a 
scale- invariant spectrum of perturbations '34*1 ) . 

For the physically reasonable models with n ^ 1 the evolution of 0, and hence of G, can be significantly altered at 
early times from what is generally assumed to be the case in scalar-tensor theories of gravity. For the case u) > —3/2, 
the value of G can be made to diverge to infinity or to zero as the initial singularity is approached, independent of 
whether or not there was an early scalar-dominated phase to the universe's history. For the case uj < —3/2, the value 
of G evolves smoothly through the bounce in the scale-factor, and is again independent of whether or not there was 
a scalar dominated phase. More complicated evolutions of <j> can also be constructed, as can be seen in Figure |31 

VIII. DISCUSSION 

We have considered spatially-flat FRW universes in scalar-tensor theories of gravity where energy is allowed to be 
exchanged between the Brans-Dicke scalar field that determines the strength of gravity and any perfect-fluid matter 
fields in the space-time. We have presented a prescription for integrating the field equations exactly for some unknown 
function A which describes the rate at which energy is exchanged. For the case of A being a linear function of r we have 
found the general solutions to the problem and for the case of A being a non-linear power law function of t we have 
been able to find a wide class of exact solutions. These solutions display behaviours that can deviate substantially 
from the corresponding solutions in the standard case, where the exchange of energy is absent. Depending upon the 
values of the parameters defining the theory and the exchange of energy, deviations in the evolution of a and 4> can 
occur at both early and late times, providing a richer phenomenology than is available in the standard theory. 

We have found that the parameter n must be bounded by the inequality n ^ 1 if the Universe is to be expanding 
at late times. For n = 1 we have found late-time power-law attractor solutions which can be used to constrain the 
parameters C2 and lu. For n < 1 we have seen that the late-time evolution will be that same as in the standard 
Brans-Dicke case, and so is subject to the same observational constraints as these theories. The parameters C3 and 
T2 have been shown to be influential only in the vicinity of the initial singularity, or at the minimum of expansion 
in non-singular solutions. These parameters are therefore less accessible to constraint by late-time observations (see, 
however, 27] where the influence on primordial nucleosynthesis is used to constrain T2). The parameter uj is, as 
always, subject to the very tight solar system constraint w > 40000 to 2(t 0]. 

These results could be of interest in attempting to explain why G is so small in the present day Universe compared 
to the proton mass scale [Gmi^ ~ 10~^^). In these models the value of G can decay away by a coupling between 
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the scalar field (j) and the matter fields which allows energy to be transferred. The small value of G is then due to 
the age of the Universe. It remains to see whether or not the late-time modifications found above are consistent with 
observations of the primordial abundance of light elements, microwave background formation and other late-time 
physical processes. Whilst being beyond the scope of this article, these studies should be able to be performed in an 
analogous way to the ones that already exist for the standard Brans-Dicke theory. 

Using the late time solutions that have been found it is possible to comment on the case of FRW cosmologies with 
non-zero spatial curvature. At early times it is expected that the effect of any spatial curvature on the evolution 
of a{t) should be negligible. From the solution H2U|I we can see that spatial curvature will dominate the late time 
evolution if the condition 

2 + 2{2 - -f)oj + 2c2 ^ 
4 + 3juj{2 - 7) - 2c2(7 - 67 - C2) 

is satisfied. If this condition is not satisfied then the power-law solution H2U|) will be an attractor as t ^ 00 even in the 
case of non-zero spatial curvature, offering a potential solution to the flatness problem. This behaviour corresponds 
to power-law inflation and it can be seen that the condition for equation (|20|l to dominate over the spatial curvature 
at late times is, indeed, also the condition that power-law inflation should occur. 

In conclusion, we have found that a direct coupling between (j) and the matter fields in scalar-tensor cosmologies 
provides a richer frame-work within which one can consider variations of G. We have shown that it is possible to 
construct models where the late-time violations of the equivalence principle can be made arbitrarily small (for A oc r) or 
are attracted to zero (for A oc r" where n < 1). This enlarged phenomenology is of interest for the consideration of the 
four-dimensional cosmologies associated with higher-dimensional theories as well as for more general considerations of 
the variation of G and its late-time value. This study has been limited to scalar-tensor theories with constant coupling 
parameters, to flat FRW cosmologies, and to special cases of cr° that allow direct integration of the field equations. 
Obvious extensions exist in which these assumptions are partially or completely relaxed. 
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